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Abstract 

Lattice models for the second-order strain-gradient models of elasticity theory 
are discussed. To combine the advantageous properties of two classes of second- 
gradient models, we suggest a new lattice model that can be considered as a discrete 
microstructural basis for gradient continuum models. It was proved that two classes 
of the second-gradient models (with positive and negative sign in front the gradient) 
can have a general lattice model as a microstructural basis. To obtain the second- 
gradient continuum models we consider a lattice model with the nearest-neighbor 
and next-nearest-neighbor interactions with two different coupling constants. The 
suggested lattice model gives unified description of the second-gradient models with 
positive and negative signs of the strain gradient terms. The sign in front the 
gradient is determined by the relation of the coupling constants of the nearest- 
neighbor and next-nearest-neighbor interactions. 

PACS: 62.20.Dc; 61.50.Ah 


1 Introduction 

Elastic deformations of materials can be described by a microscopic approach based on 
lattice equations in El and by a macroscopic approach based on continuum equations 
OT- Continuum equations for elasticity can be considered as a limit of lattice dynamics, 
where the length scales of infinitesimal continuum elements are much greater than that 
of inter-particle distances in the lattice [5]. The theory of nonlocal continuum mechanics 
was formally initiated by the papers of Eringen [6]. Nonlocal elasticity theory is based 
on the assumption that the forces between material points are a long-range type, thus 
reflecting the long-range character of inter-atomic forces [7]. Nearest-neighbor and next- 
nearest-neighbor interactions for lattice particles are important type of interactions for 
materials with non-local properties. Non-local continuum mechanics can be considered 
by two different approaches [8j: the gradient elasticity theory (weak non-locality) and the 
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integral elasticity theory (strong non-locality). This article focuses on gradient models 
of non-local elasticity suggested in [9] [HE [Hi IS] [13] - Mindlin [9] [HE H] presented a 
theory of elasticity for materials with microstructure, where it was proposed to distinguish 
quantities of the microscale and the macroscale to take into account a weak non-locality. 
Eringen [T2j also formulated a theory of nonlocal elasticity, where the integrals are replaced 
by gradients. Aifantis [13] suggested to extend the linear elastic constitutive relations with 
the Laplacian of the strain. Usually distinguish two following classes of gradient models 
with different signs of the strain gradient terms. 

The first class is formed by Laplacian-based gradient elasticity models that are de¬ 
scribed by the linear constitutive relations 
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where is the strain, a tJ is the stress, C^i is the stiffness tensor, A is the Laplace 
operator and l is the scale parameter. For l 2 = 0, we have the classical case of the 
linear elastic constitutive relations that is called the Hooke’s law. The main motivation 
for using the gradient elasticity of the form (JT]) has been the description of dispersive 
wave propagation through heterogeneous media. For the second-gradient models that are 
defined by P), it is found that the model becomes unstable for a limited number of wave 
lengths, while in dynamics, instabilities are encountered for all shorter wave lengths. The 
corresponding equation for the displacements is unstable for wave numbers k > l/l 2 . 

The second class of Laplacian-based gradient elasticity models is described by the 
linear elastic constitutive relations of the form 
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This equation has the format of equation (QE) although the sign of the higher-order term 
tends to be negative, not positive as in equation ©• Note that whereas the strain 
gradients with positive sign are destabilizing, the strain gradients with negative sign 
are stable. The strain gradients in equation ([2]) are equivalent to those derived from the 
positive-definite deformation energy density, and therefore the strain gradients in equation 
()2|) are stable. The opposite sign of the strain gradient term in equation P) makes this 
term destabilizing. Instabilities manifest themselves in dynamics by an unbounded growth 
of the response in time without external work. Instabilities are also related to loss of 
uniqueness in static boundary value problems. 

The lattice models are very important in the elasticity theory Da. At the same time, 
it is formed the opinion that the continuum models described by equation (J2]) cannot be 
obtained from lattice models. It is usually assumed that the second class of the gradient 
models does not have a direct relationship with discrete microstructure and lattice models. 
This opinion is based on the properties of the Taylor series that is used in homogenization 
procedure. This problem is described in Section 2 of this paper. In next sections, we 
propose a lattice model that allows us to remove the lack of the second class of models. 
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Moreover the suggested type of interaction allows us to have united approach to describe 
lattice models of the strain-gradient elasticity of two classes. The suggested lattice models 
give unified description of the second-gradient models with positive and negative signs of 
the strain gradient terms. A feature of suggested approach is the existence of an operation 
that transforms puna ei (see also [16], H7| ) the set of equations for coupled individual 
particles of lattice into the equation of non-local continuum. 


2 Homogenization procedure and Taylor series ap¬ 
proach 

In order to keep this paper self-contained, we briefly reproduce a derivation of the con¬ 
tinuum equation for the gradient elasticity by homogenization approach [[Til 139, 40] . 

In this section, the strain gradient models ([TJ) are derived by means of homogenization 
of the displacement field of a discrete model. In the lattice model, the particles are 
replaced by individual masses. The interactions of particle are modeled by springs that 
connected the point masses. For simplicity, it is assumed that we have one-dimensional 
lattice, where all particles have the same spring stiffness K, the particle mass M and the 
inter-particle distance d. 

The gradient elasticity models ([1]) have been derived from the continualization of the 
response of a lattice. To illustrate this approach, we will consider the one-dimensional 
system of particles and springs. All particles have mass M and all springs have spring 
stiffness K. The equation of motion of the particle n is 

MD 2 u n {t) = K ■ (u n+ 1 - 2 u n (t) + u n -i(t)) + F(n ), (3) 

where M and K are the particle mass and the spring stiffness, respectively, both of which 
are assumed to be uniform. 

In the homogenization procedure, it is assumed the continuous displacement u(x,t ) 
equals to the lattice displacement u n (t) at particle n by u n (t) = u(nd, t ), where the particle 
spacing is denoted as d. The displacement at the neighboring particles is found by means 
of a Taylor series as 

d 2 d 3 d 4 

u(x ±d,t) — u(x, t) ±d D l x u(x) + — D 2 x u(x ) ± — D 3 x u(x ) + — D x u(x ) + 0(d 5 ). (4) 

2 o 24 

Next, the displacements of the lattice medium u n ±\ (t) are expressed in terms of the 
continuous displacement. These terms are substituted into equation (J3J). After division 
by the cross-section area of the medium A and the inter-particle distance d it is found 
that 

p D 2 u(x, t) = E (d 2 x u(x, t) + ^ D x u(x, t )) + f(x) ( 5 ) 
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with the mass density p = M/Ad , the Young’s modulus E = Kd/A , and f(x) = F(x)/Ad. 
Note that all odd derivatives of u have cancelled and signs in front of second-order and 
fourth-order derivatives coincide. This equation can be written as 

d' 2 C 2 1 

D/u(x, t) = C/ D 2 x u(x, + D$u(x, t) + -fix), (6) 

where C e = \JE / p is the elastic bar velocity. 

When the kinematic relation e = D].u is used, and the equation of motion of the 
continuum is expresses as 


p D%u(x, t) = Dl<r{x, t) + f(x), (7) 

the constitutive relation can be retrieved as 

a = E (e + ^ D lej. ( 8 ) 

The positive sign in the relation (JHJ) follows directly from the positive sign of d 4 -term in 
the Taylor series (jTlh 

The second-gradient term is preceded by a positive sign, where d = Z\/l2. This pro¬ 
cedure illustrates the close relation between the discrete microstructure and the gradient 
non-local continuum with positive sign in (Jl]) and (l8j) . 

The homogenization procedure as shown above uniquely leads to a second-order strain 
gradient term that is preceded by a positive sign. The second-gradient model with negative 
sign 

a = E(e-l 2 Dle) (9) 

cannot be derived from a microstructure of lattice particles by this homogenization pro¬ 
cedure. 


3 Equations of motion for lattice particles 

In this paper, we shall use a simplest model to describe the lattice vibration, where all 
particles are displaced in one direction. We also assume that the displacement of particle 
from its equilibrium position is determined by a scalar rather than a vector. This model 
allows us to describe the main properties of a vibrating lattice by using simple equations. 

Let us consider a one-dimensional lattice system of interacting particles that are de¬ 
scribed by the equations of motion 

d 2 u (t) + °° + °° 

M—AE— = g 2 ^ K 2 (n,m ) u m (t) + g A ^ K 4 (n,m) u m (t ) + F(n), (10) 
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where M is the mass of particle, u n (t ) are displacements from the equilibrium, g 2 and <74 
are coupling constants. The terms F(n) characterize an interaction of the particles with 
the external on-site force. Let us note some properties of the coefficients K s (n,m). If 
we assume the lattice to be displaced as a whole: u n (t) = u = const, then the internal 
lattice state cannot be changed in case of absence of external forces F[n) = 0. As a result, 
equations (HH give 

+OO +OO 

£ K s (n,m)= ^2 K s (m,n) = 0, (s = 2; 4) (11) 

m=—oo m=—oo 

for all n. This requirement is well-known and it follows from the conservation of total 
momentum in the lattice [ 2 j. 

It seems that two terms of (HUD with K 2 (n, m) and K 4 (n, rn) can be combined into 
one without loss of generality. These terms are presented separately to use two different 
interaction constant g 2 and <74 for two type of interactions such as the nearest-neighbor and 
next-nearest-neighbor interactions. The advantage of such representation is manifested 
in Section 5. Using two different coupling constants allows us to derive the constitutive 
relation ([!]) and ([ 2 ]) with positive and negative depending on the relative values of the 
coupling constants r / 2 and < 74 . 

For an unbounded homogeneous lattice, due to its homogeneity the matrix K s (n,m ) 
has the form K s (n,m ) = K s (n — m) where s = 2; 4. In a simple lattice each particle is 
an inversion center, and then we have 


K s (n — m) = K s (m — n ) = K s (\n — m\). 
Using the condition (ITT]) , we can represent equations (fTOj) in the form 


M- 


dt 2 


/ , 
m=—oo 


These equations of motion take into account the translation invariance condition of a 
lattice structure with respect to its displacement as a whole. In equation (|T 2 |) the inter¬ 
action terms are translation invariant. The non-invariant interaction terms lead to the 
divergences in the continuous limit [5]. 

In this paper we consider the lattice model (fT 2 j) with the nearest-neighbor and next- 
nearest-neighbor interactions only. Then the first term /t 2 (n — m) describes the nearest- 
neighbor interaction by 

A 2 (n Tfl) 1 T $m—n, 1^ (13) 

with the coupling constant g 2 . The second term K 4 (n — m) describes the next-nearest- 
neighbor interaction by 

K 4 {n -m) = - (Sn-mp + S m _ n!2 J, (14) 
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and g 4 is the coupling constant of this type of interaction. 

We can give a discrete mass-spring system that corresponds to the suggested lattice 
model (fT2]l . In Figure 1 , we present the mass-spring system with the nearest-neighbor 
and next-nearest-neighbor interactions. 


#4 9 4 



Figure 1: Discrete mass-spring system with stiffness coefficients k 2 = g 2 and = g 4l the mass 
M and the distance d that correspond to the lattice model (fT2]i with the nearest-neighbor and 
next-nearest-neighbor interactions. 

4 Map of lattice in the continuum 

In this paper to derive a continuum equation from the lattice model we use the approach 
suggested in [hf] [15] instead of the homogenization approach described in Section 2. 

In this section, we define the map operation mmn that transforms the equations 
of motion for u n (t) of lattice model into continuum equation for a scalar held u(x,t) that 
describes displacement. 

In order to obtain a continuum equation for a lattice equation, we assume that u n {t) are 
Fourier coefficients of some function u(k,t). We define the held u(k,t) on [—K 0 /2, K 0 /2] 
by the equation 

+oo 

u{k,t)= ^ u n{t) e~ lkXn = JA{u„(f)}, (15) 

n=— oo 

where x n = nd and d = 2 it/K 0 is distance between equilibrium positions of the lattice 
particles. The inverse Fourier series transform is dehned by 

-i p+Ko/2 

u n (t) = — / dk u(k,t ) e lkXn = J 7 ^ 1 {u(k, t)}. (16) 

- tt -0 A ' q /2 
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Equations (1T5|) and (fT6l) are the basis for the Fourier transform, which is obtained by 
transforming from lattice variable to a continuum one in the limit d —> 0 (—> oo). The 
Fourier transform can be derived from (fl5|) and (fl6l) in the limit as d —)■ 0. We replace 
the lattice function 

/ \ 2 tt . . 

Un[t) = —U{X n ,t ) 

-ff o 

with continuous u{x,t ) while letting 


27m 

x n = nd = —-> x. 

Ao 


Then change (d —* 0 or A ' 0 —> oo) the sum to an integral, and equations (TT5|) and (TT6T) 
become 

/ +oo 

dx e~ lkx u(x,t) = Tluixjt)}, (17) 

-OO 

1 /»+oo 

u(x,t) — — I dk e lkx u(k,t) = J r ~ 1 {u(k 1 t)}. (18) 

2vr J -oo 

We assume that 

u(k,t ) = Limit ii(k, t), 

where Limit denotes the passage to the limit d —> 0 {K 0 —» oo). We note that u(k,t) is 
a Fourier transform of the field u(a;,t), and u(k,t) is a Fonrier series transform of u n (t), 
where we can use u n {t) = (2 tt/K 0 )u{nd, t). The function u(k,t ) can be derived from 
u(k,t) in the limit d —» 0 . 

As a result, we define the map from a lattice model into a continuum model by the 
following operation mm- 

(1) The Fourier series transform: 

•FlC\ : U n {t) ->• = u(M)- (19) 

(2) The passage to the limit d —> 0: 

Limit : u(k,t ) —» Limit{u(A;, f)} = u{k,t). (20) 

(3) The inverse Fourier transform: 

F^ 1 : u(k,t) —» 7-'~ 1 {u(fc,f)} = u(x,t). ( 21 ) 


Diagrammatically this set of operations for transformation of the displacement can be 
represented by Figure 2. 



Figure 2: Diagrams of sets of operations for displacement. 











We performed similar transformations for differential equations to map the lattice 
equation into an equation for the elastic continuum. We can represent these sets of 
transformations of the differential equations in the form of the diagrams presented by 
Figure 3. 


Equation 


From Lattice 


Equation 

for u n (t) 


to Continuum 

> 

for u(x, t) 
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Inverse Fourier 

transform T\ 


integral transform 
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Equation 


Limit 


Equation 

for u(k, t ) 


d — y 0 


for u(k,t ) 


v _y v_y v_y 


Figure 3: Diagrams of sets of operations for differential equations. 


The combination of these operations J- 1 Limit J-/\ allows us to realize a map of lattice 
models of interacting particles to models of elastic continuum. 


5 Lattice with nearest-neighbor and next-nearest-neighbor 
interactions 

Let us consider a lattice with the nearest-neighbor and next-nearest-neighbor interactions. 

The nearest-neighbor interaction is described by first terms on the right-hand side of (TT2T) 
with (fl3l) in the form 

+oo 

E K ^ n ~ m ) ( u n(f) - Um(*))= U n+1 (t) - 2 U n (t) + W n _i(t). 

m =—oo 
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( 22 ) 


















































The next-nearest-neighbor interaction in (TT2l) with (TT4l) is described by the term 


+oo 


E K t (n - m) (u n (t) - u m (t)j= u n+2 (t) - 2u n (t) + u n - 2 (t). 


( 23 ) 


Here g 2 and g A are coupling constants of the nearest-neighbor and next-nearest-neighbor 
interactions respectively. In general, we have two different coupling constants. 

The corresponding continuum equation can be obtained in the limit d —» 0 by the 
method suggested in [Tf, [T51I5]. We have the following statement regarding the suggested 
lattice model. 

Proposition. In the continuous limit (d 0) the lattice equations of motion 


M 


du r 
dt 2 


— 92 ■ ^^n+i 2 u n T Un—i'j T g 4 • |_2 2 u n T u n — 2^ T F(j 


(24) 


are transformed by the combination T 1 Limit T a of the operations PHUIill) into the 
continuum equation: 


d 2 u(x, t) 
dt 2 


d 2 u(x,t ) d 4 u(x,t ) 
dx 2 4 9a; 4 


+ 

P 


(25) 


where 


G 2 


(92 + 4gq) d~ 

M 


G a 


{.92 + 16 g A ) df 
12 M 


(26) 


are finite parameters, and f{x) = F{x)/{Ad ) is the force density, p = M/(Ad ) is f/ie 
mass density. 

Proof. To derive the equation for the field u{k,t ), we multiply equation (|24l) by 
exp {—iknd), and summing over n from — 00 to + 00 . Then 


+OO 

E 


7O +OO 

— iknd® 

e lknd ^- = g 2 ■ Y 


dt 2 


e lknd ( u n +1 — 2 u n + u n - 1 ) + 


9 4 • 


+00 

E 


+OO 

e lknd {u n+ 2 — 2w n + u n — 2 ) + 


—iknd 


F(n). 


(27) 


The first and second terms on the right-hand side of (1 T7\ are 

+OO 

lL n +\ 2 U n + U n —1 j “1“ 9 4 * ^ ^ ^ ( ^n+2 2li n + U n —2 


X) 6 

n=—oo 


+00 


+ OO 


+ OO 


= 92-1 E '= '‘”""-11 


2 E 


e~ iknd u 7 


+ Y, e- iknd u n -1 I + 
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+oo 


+oo 


+oo 


+94 • V e- lknd u n+2 


+oo 

= 92- +“ J2 

m =—oo 


e-' kmd u m 


2 E 

n=—oo 

+oo 

-2 E 


e~ iknd Ur 


e~ ikn \ 


e~ iknd u n . 2 | = 


+oo 


+oo 

= 94-1 e 2 " v 

m =—oo 


+oo 


e~ ikmd u r 


2 E 


—iknd n . i ^—2 ikd 


e~^' M u n + e 


+ e -iM ^2 

j=- OO 
+oo 

E 


I + 


-ikjd . 


e u 


j i • 


(28) 


j=-oo 


Using the definition of u(k,t), equation (l28|i gives 

g 2 . f e lkd u(k , f) — 2 u(k, t ) + e~ lkd u(k, t )) + 


+34 • (e 2ifcd u(M) - 2u(fc,t) + e~ 2 ikd u(k,t)j = 

= 92 ■ (e ikd + e~ ikd - 2)u(fc, f) + ,g 4 ■ (e 2ifcd + e- 2iM - 2)«(fc, t) = 

= 2 ^32 • (cos ( kd ) — 1 ) + 34 • (cos (2 kd) — 1 ) j h(fc, f) = 

= 2 ^-232 • sin 2 (^pj - 834 • ^sin 2 ^y^) - sin 4 (y))) ^ k,t ) = 

= - 4(32 + 434) • sin 2 ^y^) «(M) + !6 34 • sin 4 ^y^) «(M)- 
Substitution of (l29li into (1271) gives 


(29) 


M — = -4 (32+434)-sin 2 (y] “(M) + 16 34-sin 4 P] u(k,t)+J r A {F(7i)}. (30) 


2 J 


Using the asymptotic behavior of the sine in the form 


. ,'kd\ kd 1 / kd\ 3 r lN 

sll 4 t) = t _ 5 (y) +0((w)) - 


(31) 


we have 


sm 2 ,+) = Mi _ 2 1 “M! + 0 ((w) 5 ) = hf -- ++ + o( ( w)-), (32) 


W(+)=h+ + 0 ((H«. 


(33) 
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As a result, we can use the representation 


kd\ 


4 (.92 + 4.94) • sin 2 —- + 16 g 4 • sin 4 — = 


2 J 


kd\ 


2 J 


t s (kd ) 2 , , „ , (fcrf ) 4 ( kd ) 4 

= -4 (92 + 494 )- 7 - h 4 (92 + 494 ) • —— + 16 94 ■ 


48 


= -(92 + 494) • (kd) 2 + — (92 + I694) • (kd) 4 . 
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( 34 ) 


Using the finite parameter C 2 = K d 2 /M, the transition to the limit d —> 0 in equation 
gives 


where 


P = 


M 

~Ad y 


d 2 u(k, t) 
dt 2 

Kd 


= —G2k 2 u(k, t) + G^k 4 u(k, t) + J r {F(x)}, 


E = 


A ’ 


Go = - = 


E (92 + 494) d 2 


The inverse Fourier transform T 1 of 


p M 

has the form 


Ga = 


(92 + I694) d 4 


12M 


(35) 


(36) 


8t 2 


= -G 2 E~ L {k 2 u(k, t)} + G 4 E~ L {k 4 u(k, t)} + -f (x). 


Then we can use the relation T 1 {u(k,t)} = u(x,t) and the connection between the 
derivatives and its Fourier transforms 


— 1 r 7 2 ~/ 7 d 2 u(x,t) 4. 4 _. d 4 u(x,t) 


E- l {k 2 u(k,t)} = 


dx 2 


J 7 {k u(k,t)} — 4 


As a result, we obtain the continuum equation 


dx A 

This ends the proof. 


(37) 


The correspondent principle and relations (T36|) gives 

92 + 494 = K. 


(38) 


It is easy to see that the sign in front the gradient is determined by the value of the 
coupling constant 94 for next-nearest-neighbor interaction. If we have the inequalities 
94 < 0 and I 94 1 > (1/12) K, then the sign of the constant G 4 will be negative. 

If we use the kinematic relation e(x,t) = du(x,t)/dx, and the equation of motion of 
the continuum in the form 


d 2 u(x,t) lda(x,t) 1 , 

—dF~ = p-aT- + p J {x) " 


(39) 
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then the constitutive relation is represented as 


a(x, t) = E 


( (92 + 4 g 4 ) p d 2 (g 2 + lQg 4 )pd 3 4 d 2 e(x,t) \ 

v EM £[X ’ ’ 12 EM dx 2 ) 


with the mass density p = M/Ad, and the Young’s modulus E = Kd/A. Using 
rewrite relation (UOll in the from 


a(x, t) = E 


(e{x,t) + 


(.K + 12g 4 ) d 2 d 2 e{;x, t ) \ 
Y2K dx 2 ) ' 


(40) 


(41) 


The second-gradient term is preceded by the sign that is dehned by sgn(/i + 12g 4 ). The 
scale parameter l 2 of the gradient elasticity is connected with the coupling constants of 
the lattice by the equation 


K + 12 (/ 4 | d 2 
12 K 


(42) 


Equations (j35|) and (jH]) give the close relation between the discrete microstructure of 
lattice and the gradient non-local continuum. 

The sign in front the gradient in d4T]) is determined by the relation of the coupling 
constants for nearest-neighbor and next-nearest-neighbor interactions. We can list all 
possible cases: 

(1) If the coupling constants of the next-nearest-neighbor interaction is 


1 1 

— 7 92 < 9a < —— 92, 
4 lo 


(43) 


l2 8 2 £(x,t) 
dx 2 


then we have the stress-strain constitutive relation with nagetive sign 

a(x, t) = E ^e(a;, t) — l 2 
(2) If we have the inequality 

1 

94 > —— < 72 , 

16 

then the sign in the constitutive relation is positive 

a(x, t) = E ( e(x, t ) + l 


2 d 2 e(x,t ) 


dx 2 


(44) 


(45) 


(46) 


(3) If we have g 4 < —(1/4) g- 2 , then the constant G 2 will be negative and the sign of 
G 4 is the same as one of g 4 . 

(4) If we have c/ 4 = —(1/4) g- 2 , then the constant G 2 and K are equal to zero and the 
sign of G 4 is the same as one of g 4 . 
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(5) If we have the condition 


9a = ~Yq 92 ' ( 47 ) 

then the constant G 4 is equal to zero and we have the well-known Hooke’s law without 
weak nonlocality 

a(x,t) — E e(x,t). (48) 

As a result we can state that the suggested lattice model gives us an unified approach 
to the second-gradient elastic continuum models with positive and negative signs of the 
strain gradient terms. The close relation between the discrete microstructure of lattice 
and the gradient non-local continuum is given by equations (j38|) and (l42]h 

The proposed lattice model with the nearest-neighbor and next-nearest-neighbor in¬ 
teractions uniquely leads to second-order strain gradient terms that are preceded by the 
positive and negative signs. The lattice models with positive value of coupling constant 
$4 of lattice model leads to the continuum equation with the positive sign in front of 
the parameter l 2 . This continuum equation is unstable for wave numbers k > l/l 2 . The 
instability leads to an unbounded growth of the response in time without external work. 
The negative value of coupling constant g 4 of lattice model can lead to stiffness coefficient 
of the next-nearest-neighbor interaction with non-convex elastic energy potentials in the 
discrete mass-spring system. At the same time, the correspondent continuum equation 
with the negative sign in front of the strain gradient are equivalent to those derived from 
the positive-definite deformation energy density, and therefore these continuum models 
are stable. We should note that continuum limits for discrete and lattice systems can 
be considered without convexity hypotheses on discrete energy densities [18]. In addition 
there exist metamaterials with negative stiffness. In the plastic deformation range one 
can observe a decreasing part of the displacement-force curve, where stiffness along this 
part of the curve is negative. Viscoelastic materials, nanofilms and molecular chains con¬ 
taining a negative-stiffness phase, anomalies in stiffness and damping have been observed 
experimentally HSHUUZEJI 22 ]. One of the negative stiffness sources can be obtained from 
phase transforming materials in the vicinity of their phase transition. A theoretical de¬ 
scription of the underlying mechanism from a microscopic viewpoint has been suggested 
in [23]. Drugan [24] demonstrated that elastic composite materials having a negative 
stiffness phase can be stable. A stability analysis for elastic composites with non-positive- 
definite phase was made by Kochmann and Drugan in [2S]. Lee and Goverdovskiy [26] 
demonstrated that negative-stiffness elements can be used to improve vibration isolation 
systems. Elastic composite materials with negative stiffness inclusions can have extreme 
overall stiffness and mechanical damping [ 2 T 1 f 28 ] . 
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6 Conclusion 


In this paper lattice model with the nearest-neighbor and next-nearest-neighbor interac¬ 
tions for strain-gradient elasticity of continuum is suggested. This model can be considered 
as a microstructural basis of description for two types of gradient models with positive 
and negative signs. 

We can formulate the main results of this paper in the following short form: (a) The 
method, suggested in [T4J US] and using the Fourier series and integral transformations, is 
expanded to elasticity theory; (b) This method is applied to modified lattice model with 
two coupling constants; (c) The gradient elasticity models with negative sign is derived 
from this lattice model that have not be done before. Let us give some details about the 
suggested lattice model with two coupling constants and the proposed method that uses 
the Fourier transformations 

In the lattice models of gradient elasticity we suggest to use two different coupling 
constants g 2 and g±. It allows us to derive the constitutive relation that contains these 
constants. The usual method that based on the direct Taylor expansions (0]) allows us 
to derive the constitutive relation (jS]) with positive sign only. This traditional approach 
corresponds to the case g± = 0, where the next-nearest-neighbor interactions are not taken 
into account. It cannot give the elastic constitutive relations (J2J) with minus sign. Only 
incorporation of two types of interactions with two different constants for the nearest- 
neighbor and next-nearest-neighbor interactions may result in the continuum limit to the 
gradient elasticity theory with positive and negative signs. 

The method based on the Fourier series and integral transformations is represented 
by Figure 3. This method has been suggested in the papers mm , and it was applied 
in m m- Unfortunately this method is not applied to models of gradient elasticity. 
Moreover there is a common opinion that the gradient models with a negative sign in the 
constitutive relation ((2]) cannot be obtained from lattice models. In this paper, the Fourier 
method is expanded to lattice model with two coupling constants, which is represented 
by Figure 1. It allows us to prove that the gradient continuum model with negative sign 
can be obtained from microscopic (lattice) model. This is one of the main results of this 
paper. We also suggest the conditions (1T3D for coupling constants, when the gradient 
continuum model with negative sign can be obtained by the lattice approach. 

Let us note some possible extensions of the suggested lattice approach. Using the 
second-, third, ... nearest-neighbor interactions in the lattice in addition to the nearest- 
neighbor interaction it can be easily generalized for the case of the high-order gradient 
elasticity. For this generalization, we can use the terms K 2j (n, rn) = — + 

where j = 1,2,3,..., that describe an interaction of the n-particle with two particles 
with numbers n ± j. Using models of lattice with long-range interactions, we can get 
continuum equations with the fractional derivatives [29j of non-integer orders by the 
methods suggested in [14], [15]. Therefore the suggested lattice model can be extended 
to get equations for elastic continuum with power-law non-locality [30, SE [32, |M1 33]. 
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We also assume that the proposed lattice approach can be generalized for lattice with 
fractal dispersion law (351 ESI E7]. ft allows us to get a possible microscopic basis for 
gradient elasticity of fractal materials [38] that can be described by different tools (see 
for example [3H HZ]) - If we consider lattice equation (1TQ|) . (1T2|) . (1241) with nonlinear site 
force F(u n (t)) instead of F(n), then we can get differential equations with nonlinear term 
f{u(x,t)) for a continuum in the continuous limit. It is well-known that the localized 
soliton solutions for lattice models are very important for nonlinear theories [131 H]. 
Solitons, its interactions, and correspondent discrete lattice models can be considered 
for different type of physical systems. For example, equations for an electromagnetic line 
with the second-order couplings can be used to describe the bound solitons with oscillating 
tails in this line (43]. The approach suggested in this paper can be generalized in order to 
consider continuum models that correspond to discrete model of a dynamical lattice with 
the on-site nonlinearity and both nearest-neighbor and next-nearest-neighbor interactions 
between lattice sites (Ml S3 HE]. We also assume that the unbounded lattice models 
suggested in this paper can be extended to describe bounded lattices and correspondent 
continuum models analogous to electromagnetic line described in |J5]. We also can assume 
that suggested approach can be useful for nonlinear deformable-body dynamics [45] for 
the case of weak non-locality. 
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